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Abstract 

Wc study scattering coefficients in black hole spacctimes using analytic properties of complexified 
Q-c wave equations. For a concrete example, we analyze the singularities of the Tcukolsky equation 

r~| ■ and relate the corresponding monodromies to scattering data. These techniques, valid in full 

generality provide insights into complex-analytic properties of greybody factors and quasinormal 

modes. This leads to new perturbativc and numerical methods which are in good agreement 

with previous results. 
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1 Introduction 

In the study of quantum field theory in curved spacetime, perhaps the most basic object worth 
considering is the two-point function in a black hole background. This quantity describes the 
amplitude for a particle to scatter off of the black hole and contains a wealth of information of 
formal and observational interest. In particular, it appears in the formula for Hawking radiation, 
thus playing a role in the description of quantum mechanical black holes; indeed, the low frequency 



behavior of scattering coefficients can be matched precisely with proposed conformal field theory 
descriptions of black hole microstates. Observationally, the poles in the scattering coefficients are 
quasinormal modes, which characterize the late-time relaxation towards a black hole equilibrium 
state. Formal properties of these quasinormal modes are related to the stability of black hole 
spacetimes, which is in many cases an open and difficult problem. 

In computing these scattering amplitudes the primary obstacle is technical — the field equations 
cannot be solved analytically. In this paper we describe techniques, based on global analytic 
properties of probes of black hole geometries, that allow us to extract certain precise features of 
scattering amplitudes even without a direct analytic solution. This approach has the advantage 
of distinguishing global versus local data, allowing us to parameterize the dependence of certain 
physical observables on analytic properties of the solution. 

Some of the results we present here are not new; we refer the reader to the reviews of [H [21 [3] 
and citations therein. Global aspects of the relevant class of ODEs have been studied to some 
extent in prior literature. See e.g. [U [6] for an extensive analysis of properties of the solutions 
to confluent Heun equations and [7J [H] for an interesting related derivation of greybody factors for 
the Kerr black hole. 

Our analysis starts with a remarkable fact: in a broad class of black hole backgrounds, including 
Kerr, the field equation separates; radial dependence is determined by an ordinary differential equa- 
tion with regular singular points at the event horizons and an irregular singular point at asymptotic 
infinity. This means that if we complexify the radial coordinate, the field is locally a holomorphic 
function of the radius with branch points at the horizons and infinity. These branch points reflect 
the fact that two linearly independent solutions of the wave equation will mix when they are trans- 
ported around a non-trivial loop in the complex radial plane. The failure of the field to be globally 
holomorphic is encoded in monodromy matrices associated with the singular points at the horizons 
and infinity. At the horizons, these matrices have a simple physical interpretation: their eigenvec- 
tors correspond to the solutions that are purely ingoing/outgoing. At infinity, ingoing/outgoing 
boundary conditions instead correspond to eigenvectors of a formal monodromy matrix (as we will 
discuss in section 12. 4p . The scattering amplitude is then encoded in the change-of-basis matrix 
between two sets of monodromy eigenvectors. The advantage of this procedure is that monodromy 
data can, in many cases, be extracted from purely local features of the geometry, without the need 
to evolve wave functions from the horizon to asymptotic infinity. The topology of the complex 
radial plane (with singular points removed) leads to a non-trivial global constraint, allowing us to 
extract certain features of scattering data without ever solving a differential equation. 

Although elegant, this algebraic procedure does not completely fix the reflection and transmis- 
sion coefficients. In simple backgrounds, such as BTZ black holes, it allows us to compute the norm 
of the transmission and reflection coefficients — i.e., the greybody factor — but not their phases. The 
computation of quasinormal modes in these simple cases, then, only requires a little more work. 
For more complicated backgrounds, such as Schwarzschild or Kerr black holes, the singular point 
associated with asymptotic infinity is irregular, so it exhibits Stokes phenomena; this allows global 
features of the black hole geometry to sneak into monodromy data at infinity. Describing these 



data thus requires either perturbative or numerical computations, which we develop. Beyond this 
complication, as will become clear, there is another subtlety that prevents us from precisely com- 
puting scattering data in this manner. Nevertheless, our method provides new analytic insight. 
For example, quasinormal frequencies are expressed as roots of a certain transcendental equation, 
allowing us to qualitatively understand their analyticity properties in the complex frequency plane. 
In the following section, we describe the monodromy technique, focusing on the extraction of 
scattering data and a treatment of Stokes phenomena. In section El we apply this to the Kerr black 
hole. In section U we consider the symmetries of the relevant wave equation which, when combined 
with monodromy techniques, allow us to develop new methods for the study of Kerr quasinormal 
modes. Many of the technical results are relegated to a series of appendices. We also refer the 
reader to a companion paper [S], which relates this technique to microscopic treatments of black 
hole quantum mechanics in terms of a conformal field theory. 

2 Overview of the Monodromy Technique 

Our goal is to study the dependence of scattering coefficients and quasinormal modes on analytic 
(global) properties of linearized fluctuations around a black hole background in asymptotically flat 
spacetime (other asymptotics can be treated in the same vein, with appropriate modifications). A 
typical scattering computation requires one to solve a wave equation with specified boundary con- 
ditions. The essential point is that when those boundary conditions are specified at singular points 
of the wave equation, then they will be intimately connected to analytic properties of solutions 
when the independent variable is complexified (typically the radial coordinate). In this section, 
we will explain how knowledge of the analytic properties of the solutions can be used to compute 
scattering coefficients, reducing the computation to a simple exercise in linear algebra. 

This section is rather technical; an eager reader more interested in black hole physics is invited 
to skim this section and head directly to section [3j A more complete treatment of the relevant 
ODEs can be found in, e.g., ^OlfTT] . 

2.1 Preliminaries: ODEs and Monodromies 

We will focus on second order ordinary differential equations of the form 

8 Z (U(z)d z il>(zj) - V(z) if>(z) = . (2.1) 

Such equations can always be rewritten as two coupled first order ODEs: for instance, introducing 
a two-component column vector ^ satisfying the ODE 





^2 



=: A{z)V , (2.2) 



which reproduces (|2.ip when we identify \&i with ift since then ^2 equals U{z)d z ip. The space of 
solutions to this ODE is two-dimensional, so we can choose a linearly independent basis of solutions 
vl/W and \I / ( 2 ) and collect them into a fundamental matrix 

$( z ) := ( q,W ^(2) \ . (2.3) 

The linear independence of the two solutions is equivalent to the invertibility of &(z). 
The Wronskian of two solutions, ipi(z) and ip2(z), is 

W(ihM ■= U{z)faz)d x M*)-Mz)dzih(zJ) • (2-4) 

Note that the determinant of &(z), which equals the Wronskian of ^f\ with ^\ , is a constant 
since tr(A) = 0. This follows from the path-ordered representation, <&(z) = V\exp(f z A)\&q, 
which shows that det($(z)) = det($o) is independent of z since 

det (v{ef ZA }) = det (e^ A ) = e^^ = 1 , (2.5) 

where V denotes path ordering. 

We now analytically continue to the complex z-plane and consider cases where A{z) is mero- 
morphic; all cases considered here will certainly satisfy this requirement. In this case, one may 
think of A(z) as a flat SL(2; C) connection^] where allowed gauge transformations must themselves 
be meromorphic functions of z. The gauge-invariant information of A(z), namely its holonomies, 
is thus encoded in the analytic structure of <&(z); in particular, if <&(z) is meromorphic, &(z) is the 
gauge transformation that sets A(z) to zero. To more directly see this relationship, follow &(z) 
around a closed loop 7 in the complex z-plane, calling the result $ 7 (ar). Since A(z) is meromorphic, 
the differential operator d z — A returns to itself, which implies that 3> 7 (z) must again be a funda- 
mental matrix for the ODE, however it need not be equal to <&(z): given one fundamental matrix, 
we can always multiply it from the right by a constant invertible matrix to obtain another (i.e., we 
can choose a different linearly independent basis of solutions). By the definition of a fundamental 
matrix, then, 3> 7 (z) must equal 4>(z)M 7 for some invertible matrix of constants M 7 : 

$ 7 (z) = V{e^ A }^{z) =: ^{z)M 1 . (2.6) 

By definition, if our loop 7 does not cross any branch cuts of <&(z) then M 7 = 1, so M 7 is a measure 
of the lack of meromorphicity of &(z) and is called a monodromy matrix. A rearrangement of (12. 6p 



M 7 = ${z)- x V{(fi A }$(z) (2.7) 

emphasizes the relationship between the holonomy of A{z) around 7 and the monodromy matrix 



Given a meromorphic A(z), not necessarily of the form in H2.2|) . one might expect it would be a GL(2; C) 
connection, however one can always perform a meromorphic gauge transformation to force ti(A(z)) — 0, as satisfied 
by (EU). 



M 7 : they are in the same conjugacy class. 

In particular, if one can find a gauge in which A has no poles enclosed by 7, then M 7 = 1. In 
this way, poles of A that cannot be removed by gauge transformations correspond to branch points 
of &(z) and are associated with non-trivial monodromy matrices. Since the conjugacy class of the 
holonomy "Pjexp (<f> A)\ is independent of variations of the loop 7 that do not cross other branch 
points of $, the monodromy matrices associated with the branch points define an embedding of the 
first homotopy group IIi(P 1 \{branch pts}) into SX(2;C), where we have compactified the z-plane 
to a P 1 by adding the point at infinity (which frequently is itself a branch point, as determined 
unambiguously by the connection A). 

This has an implication that is key to the rest of our study. Let z = Z{, for i = 1, . . . ,n, be 
the locations of all branch points of &(z), and let Mj be the monodromy matrix associated with a 
loop that encloses only the branch point at Z{. If we follow &(z) around a path enclosing all branch 
points, the other side of the loop encloses no branch points and so the monodromy around that 
loop must be trivial. In other words, 

MiM 2 ■ ■ ■ M n = 1 . (2.8) 

This innocent relation is actually quite interesting. The conjugacy class of each individual Mj 
can often be computed quite simply from local information of the differential equation (with an 
important caveat in section 12741) . while (I2.8P is a relation among these local data — it is a piece of 
global information. Computing scattering coefficients is an example of a problem where we require 
global information about our solutions (relating boundary conditions at different points in the z 
plane), so the rest of our work will explore how we may exploit this relation. 

As we alluded to above, if $ has a branch point at Zi, then A has a pole there. The converse is 
not true because it may be possible to remove the pole in A by a gauge transformation. Fortunately, 
there exists a simple algorithm for choosing gauge transformations to reduce the order of the pole 
to some minimal integral value called the Poincare rank, Rj G No, so there is a gauge where A{z) 
takes the form 

(z - z t )-^- 1 A (z) , (2.9) 

where Aq(z) has a convergent Taylor series expansion in some neighborhood containing z = Zi. A 
simple pole Rj = is called a regular singular point, while a higher order pole Rj > is called an 
irregular singular point of rank Rj . 

The distinction between regular and irregular singular points might seem artificial, but their 
implications for <I> are starkly different: regular singular points correspond to algebraic or logarith- 
mic branch points in <&, while irregular singular points correspond to essential singularities in $ 
and exhibit Stokes phenomenon. Understanding regular singular points will suffice to illustrate the 
essence of our approach to scattering computations, so we focus first on this case in sections | 
and 12. 3| discussing the added complications of irregular singular points in section 12.41 



2.2 Monodromies and Regular Singular Points 

Our goal now is to find monodromies and the local behavior of $ around singular points. To find 
the conjugacy class of Mi, first perform a gauge transformation to bring A to the minimal form 
(|2,9p . where now Rj = 0. For loops enclosing a single pole, this renders the path ordering trivial in 
the limit that the loop approaches the pole, but there can be one more subtlety: if the eigenvalues 
of Res Zi (A) differ by nonzero integers, one must perform additional gauge transformations to make 
the eigenvalues equal (for a systematic procedure, see, e.g., [TO]). Once that is done, the conjugacy 
class of the monodromy matrix (|2.7[) is easily determined to be 

Mi <* exp I - z i y 1 A {z) = e 27ciAo ^ , (2.10) 

where '— ' denotes equal up to conjugation by an element of SL(2, C). It is then a general theorem 
that around z = z; we can write 



$(*) = x>-o n $n (*-*)"*> ( 2 - n ) 



\n=0 / 

where the $ n are constant matrices with <3?o invertible, and the series has a nonzero radius of 
convergence around z%. Following the solution around Zi via (z — z{) — > e 2m (z — zi), we see that the 
monodromy matrix is 

M t = e 2wiN > . (2.12) 

The conjugacy class of Ni is thus the same as that of Ao(zi), so this can be easily read off the ODE 
(still assuming that A has been put into minimal form around Zj). 

2.3 Monodromies, Boundary Conditions, and Scattering 

The relationship between monodromies and choices of boundary conditions at the regular singular 
points — which is certainly a natural place to set boundary conditions for black holes where the 
horizon is a regular singular point — follows readily from (|2.1ip . For simplicity, suppose that Mi 
has distinct eigenvalues e Zf2wc " i (so Ni has eigenvalues ±iaj, up to shifts by integers), then we are 
free to choose our fundamental matrix to diagonalize Mi, in which case 

«.)-{* + Ol,-«))(l-f (z _°,_ m ). ("3, 

Approaching Z{ from a direction where Im(iajln(z — Zij) 7^ 0, we see explicitly that one column 
corresponds to ingoing waves and the other to outgoing waves. Diagonalizing the monodromy 
matrix Mj therefore corresponds to choosing a basis with definite boundary conditions at ZjO 



2 When Mi has a nontrivial Jordan block (hence, equal eigenvalues), it seems more difficult to identify diagonal- 
ization of Mi with ingoing and outgoing boundary conditions since one solution always has a logarithmic branch cut. 
However, depending on the particular problem being studied, other choices of boundary conditions may be perfectly 



A scattering computation often involves finding the change of basis between solutions that 
are ingoing/outgoing at one singular point and solutions that are ingoing/outgoing at another 
singular point (or, sometimes, normalizable/non-normalizable) . To be more concrete, focus on two 
singular points, z\ and Z2, and two fundamental matrices <3?i and <1>2 that diagonalize M\ and M2, 
respectively. Since there are only two linearly independent solutions to the ODE, <£^~ $1 must be 
a constant matrix and, in fact, 

AW = $^$1 (2.14) 

is nothing more than a change of basis from left-eigenvectors of M2 to left-eigenvectors of M\. So 
far we have said nothing of the normalization of the various bases of solutions, so we are free to 
rescale the columns of $1 and <I>2) meaning the change of basis .M2-S.1 is thus far defined only up 
to multiplication by a diagonal matrix on either side 

AW ~ (* J-M»-i(* J • ( 2 - 15 ) 




If, for example, we have a signature (1, 1) inner product for our solutions — the Klein-Gordon inner 
product in the case we will study — then we can choose the columns of $1 and of $2 to be normalized, 
lifting (most of) the rescaling ambiguity by forcing Af2_»i G SU(1, 1): 



AW = \ r T \ |ft|2 + \ T \2 = ! } ( 2.i 6) 



where 1Z and T are the reflection and transmission coefficients, respectively. The only remaining 
ambiguity will be in changing the normalizations by phases, but these will fortunately cancel when 
computing \TZ\ or \T\. 

We have now seen how monodromy matrices relate to boundary conditions at regular singular 
points and we have seen how to compute the conjugacy class of the corresponding monodromy 
matrices. If we knew the actual form of the monodromy matrices rather than just their conjugacy 
classes, computing the scattering matrix would reduce to the simple linear algebra problem of 
finding the change of basis between the left-eigenvectors of M2 and the left-eigenvectors of Mi, 
then using the relation (|2.15|) to fix A^^i to be an SU(1, 1) matrix. In fact, there are many 
physically interesting situations where knowing the conjugacy classes and the relation f|2.8j) allows 
one to reconstruct the matrices in a common basis. For example, when there are two singular 
points then the monodromy matrices are inverses of each other and the scattering matrix is the 
identity. However the nature of the singular point is important, so don't be misled here: for example, 
the Coulomb potential has two singular points, yet scattering is nontrivial because infinity is an 
irregular singular point and the choice of plane- wave boundary conditions at the irregular singular 
point does not diagonalize the associated monodromy matrix, as discussed in section 12.41 

A more interesting situation, relevant to our discussions of black holes, is when there are three 
singular points. In this case, knowledge of the conjugacy class of Mi, M2, and M3, together with 



sensible, for example demanding regularity of the solution at z%. 



the global relation in (|2.8p . is enough to reconstruct the matrices themselves in a common basis. 
Explicitly, given 

det(Mj) = l, tr(M 4 ) = 2cosh(2vrai), M* ^ 1 , for i = 1,2,3, (2.17) 

and equation (|2.8p . 

M 1 M 2 M 3 = 1 , (2.18) 

then a common basis is given by 

U -1 \ /2co S h(2 TO2 ) e— \ 

y 1 2cosh(2vrai) J \ -e" 2 ™ 3 J 

( e 2 ™ 3 \ 

M 3 = , . so- (2-19) 

^ 2(e- 27ra3 cosh(2vrai)-cosh(27ra 2 )) e" 2 ™ 3 y 

We can readily write down a change of basis from left-eigenvectors of M 2 , diagonalizing M 2 as 
diag{e~ 27rQ2 ,e 27ra2 }, to left-eigenvectors of Mi, diagonalizing M\ as diag{e~ 2?rai , e 27rai }, as 

/u ~ / sinhvT(a 2 -ai + a 3 ) sinhyT(a 2 + «i + a 3 ) \ 

y sinh7r(a2 + ai — a 3 ) sinh7r(a2 — «i — a 3 ) / 

Whether this can be made an SU(1, 1) matrix depends on details of the ckj, however when it is 
possible — for instance, when the ckj are all real and A^2->-i is invertible — we can read off the norm 
of the transmission coefficient without even computing the required diagonal transformation (I2.15P : 

\T\ 2 = 1 - \K\ 2 = sinh(2vTai) sinh(27ra 2 ) 

sinhvT(a 3 + ol\ — a 2 ) sinh7r(a 3 — a\ + a 2 ) 

Of course, ODEs with three regular singular points have hypergeometric functions as solutions 
and are therefore well understood, and we have verified that these formulas are correct. The 
challenge for applying the same methods to scattering off black holes becomes apparent in the 
next section where we discuss the consequences of an irregular singular point, which black hole 
backgrounds have at asymptotic infinity (a consequence of plane waves having essential singularities 
at infinity). The basic idea for computing the scattering coefficients will be the same, but there 
will be additional steps and subtleties. 

2.4 Formal Monodromies are Fake: Irregular Singular Points 

Asymptotically flat black holes have irregular singular points at infinity. This means that &{z) has 
an essential singularity at z = oo, reflecting the fact that the solutions asymptote to plane waves. 
Although the black holes we will study in later sections have singularities of Poincare rank 1, 
in this section we will consider general Poincare rank R > 1. For simplicity, let A(z) have a rank 



R > 1 singularity at z = oo, which means there is a gauge in which A(z) can be written as 

z^Aoiz) , (2.22) 

where Aq(z) has a convergent Taylor series expansion in z~ l around z = oo and ^4o(°°) 7^ 
(furthermore, there is no gauge in which A is less singular at z = oo). For simplicity, let Ao(oo) 
have maximal rank|f| and choose the gauge so that Aq(oo) is diagonal. Then there exists a formal 
fundamental matrix of the form 

$ f {z) = P(z)e A ^ , (2.23) 

where P(z) is a non-negative power series in z~ l (generally not convergent) and 

R 

A(z) = ^A a z a -A log(z), (2.24) 

a=l 

where the A a and Ao are all constant diagonal matrices, determined by requiring $j solve the 
ODE order by order in z~ l — in particular, A R = ^4o(°°)- It would seem that we can read off the 
monodromy from these solutions to be e 2n °, as before (around infinity, z — > e~ 2m z is the positive 
direction); however, this is not so because the calculation is complicated by the fact that P{z) 
is just an asymptotic series, not a convergent series. For this reason, e 27r ° is called the formal 
monodromy, but we will shortly see its relation to the true monodromy. 

For describing boundary conditions at z = oo, though, we care about diagonalizing k(z) (as 
opposed to the true monodromy) since this describes the dominant behavior of the solutions at 
z = oo, telling us whether the solutions are ingoing or outgoing. On the other hand, the quantity 
entering the product relation (j2.8[) is the true monodromy, so it is crucial to understand the rela- 
tionship between the two. The distinction and relation between them is intertwined with the fact 
that solutions of ODEs around irregular singular points exhibit Stokes phenomenon. The defining 
feature of Stokes phenomenon is that it arises when one attempts to describe one function (e.g., 
a fundamental matrix of true solutions) in terms of a function with a different branching struc- 
ture (e.g., a formal fundamental matrix). Again, this arises here because the formal solutions are 
generally not convergent series. 

To give the general idea, let us be a bit schematic, clarifying the details in the following para- 
graphs. Consider the following formal expression: 

Qfiz)- 1 ®^ = e- A(z) P{zy 1 P{z)e A(z K (2.25) 

The important observation to make is that we can only say that 

P(zy 1 P(z) ~ 1 (as z -> oo) , (2.26) 



3 This is called the non-resonant case, the resonant case takes a bit more work, essentially requiring one to consider 
the ODE on some multi-sheeted cover of the z-plane. 



in words, P(z) P(z) is asymptotic to 1 as z — > oo. This means it could differ from 1 by something 
with no series expansion around z = oo, e.g., e z , and this is precisely what happens. If a product 
such as (|2.25|) were between two actual fundamental matrices rather than two formal ones, we would 
expect the result to be a constant matrix. Since formal fundamental matrices are asymptotic to 
actual ones, the limit of (J2.25P as z tends to oo will be a constant matrix S, called a Stokes matrix, 
with components 

Sij = lim e - A » (z)+A ^ (z) (fe + 0(z- 1 )), (2.27) 

where the 0(z~ l ) terms are asymptotic to as z — > oo. Clearly, then, Su = 1, but what about 
off-diagonal entries? When we approach z = oo along a ray on which Re((— A R) jj + A R jj)z R ) < 0, 
the exponential forces Sij to vanish. On the other hand, when Re((— A a> a + A R jj)z R ) > 0, this 
may combine with the 0(z~ 1 ) terms to produce a finite result, implying that Sij has an upper 
or lower triangular structure. This is most evident very near Im((— A Ri jj + A R jj)z R ) = 0, where 
Re((— A Ri ii + A H jj)z R ) is maximal or minimal, as one can observe by studying the classic example of 
the Airy function. This is called a Stokes ray and represents the transition between two asymptotic 
expansions Q 

Now we can be more precise. Divide the neighborhood of the irregular singular point at z = oo 
into wedges, 

n k := { ^l < arg ( z ) + lArg(A R)22 - A R) n) < ^±^ } , keZ, (2.28) 

where < 5 <C it. The preceding paragraph illustrates that we should only think of a given asymp- 
totic expansion (e.g., one for which P(oo) = 1) as being asymptotic to a given true solution within 
a wedge, fJfc. Once we cross the next Stokes ray, the same true solution will have a different formal 
fundamental matrix (e.g., its asymptotic expansion in the next wedge will begin with P(oo) ^ 1). 
Let $(z) be a true solution such that on D,q 

${z)\ Uo ~ $/0O|n o (asz^oo, ze n ) . (2.29) 

To circle z = oo in the positive direction, we follow $(z) from f^o — > fi_i —)•...—)• ^_2r- Notice 
that along Q,q H 0_i, we have that Re(z R (— A Rj n + A R] 22)) > 0, so the first Stokes matrix is upper 
triangular. This means that the subdominant solution, the first column of <&, may be added to the 
dominant solution, the second column of <&. Thus, 

Hz)\ n _ ± ~ d> / (z)| n _ i (^o)=:$ / (z)| n _ i 5 (asz^oo, z € Sl-i) . (2.30) 

with Cn constant. 



4 Around this ray, one of the asymptotic solutions is maximally larger than (dominant) the other (subdominant), 
and the errors obtained by truncating the expansion of the dominant one become of the same order as the subdominant 
solution. This is classically phrased in terms of the coefficient of the subdominant solution changing discontinuously 
across the Stokes ray, but it has been shown that, in reality, the transition can be represented in a sharp but continuous 
fashion 1121. 
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At the next overlap, 0_i n Q~2, the roles of dominant and subdominant solutions reverse, so 
the relevant Stokes matrix will be lower triangular: 

Hz)\ Q _ 2 ~ QfWl^ic^So^Qfiz^S-rfo (asz^oo, z G tt_ 2 ) , (2.31) 

and the Stokes matrices Sk continue alternating between upper and lower triangular. Finally, 

S1_2r = ^0) so 

$(*)L , ~ $ /(*)|n , 5-2R+1---50 = $/(z)L e 2mA °5_ 2R+ i---5o (as z -> oo, z G n_ 2a ) 

(2.32) 
where e 27r ° is the formal monodromy. This provides us with the identification of the true mon- 
odromy as 

M oo = e 2mA °S. 2R+1 ---S . (2.33) 

The hard work comes in determining the C/-, commonly called Stokes multipliers. Their values are 
not solely determined by the local data of the singularity, they depend on all terms in the connection. 
Some useful references with a more complete discussion than presented here are [10t \TT\ PT3] . In 
appendix |Bl we review a simple numerical method to compute Cfc. 

3 Example: Kerr Black Hole 

The details in the following section will focus on the study of a scalar field in the Kerr black 
hole background, but the methods are readily extendable to a broad class of physically relevant 
situations. We will first revisit the wave equation of the probe with particular emphasis on exploiting 
the machinery presented in section [2j 

3.1 Wave Equation 

Kerr is the ideal playground. The solution contains the complexity needed to illustrate the method 
in the most simple manner. To setup our notation and conventions, we start by reviewing aspects 
of the geometry. In Boyer-Lindquist coordinates, we have 



A , „ . 9„ , ,\2 . „ ,„9 . sin 2 fl f , , . 9n , , „,2 

A"' S 



ds 2 = -dr 2 --(dt- a sin 2 Odcp) +Zd6 2 + ^^ ((r 2 + a 2 ) dcf) - adty , (3.1) 

where 



A = r 2 + a 2 - 2Mr , S = r 2 + a 2 cos 2 9 . (3.2) 

This describes a generic 4D black hole with mass M and angular momentum J = Ma. The inner 
and outer horizons are located at 

r±=M± VM 2 - a 2 . (3.3) 
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The Klein-Gordon equation for a massless scalar is 

1 



-A, 



-fj 



gg^d^) = 



(3.4) 



Generalizations to massive probes and tensor fields is straightforward. Expanding in eigenmodes 



iP(t,r,6,<l>) =e- i " t+m<t) R{r)S{6) , 
and using (|3.ip . makes equation (|3.4|) separable. The spheroidal function S(8) satisfies 



(3.5) 



sm( 



-d e (sin 6d e ) 



m 



2 2 2 

+ a J w J cos" 



8111 



S(6) = -K e S(9) 



(3.6) 



for eigenvalue Kg, and the radial equation for R(r) is then 



d r Ad r + 



(2Mr + w — am) 2 (2Mr_uj — amy 



(r — r + ){r + — r_) (r — r^)(r + — r_) 



+ (r 2 + 2M(r + 2M))w 2 



fi(r) = K,i?(r) .(3.7) 



Equations (|3.6p and (J3TTJ) are of the confluent Heun type — each with two regular singular points 
and a rank-1 irregular singularity — coupled by a separation constant Kg. 

Most of our discussion deals with the details of the radial equation. Our interest in the angular 
equation (|3.6p is solely to determine the eigenvalues Kg. This can be done systematically by treating 
the w-dependent term as a (not necessarily small) perturbation [14|. I15j . A basis of solutions is given 
by 

S(9) = Y,dTA^)Y £ T(6), m = -£,... ,£, £ = 0,1,2,..., (3.8) 



where Y^ n {6) are associated Legendre polynomials of the first kind. The coefficients d^,(aco) are 
determined by a recursive relation obtained by replacing (|3.8p in (|3.6p . A similar recursive method 
also determines Kg to be 

oo 

K e (auj) = ^2k n (auj) 2n , 

n=0 

where the first few terms are 

2m 2 + 1 - 2£{£ + 1) 



(3.9) 



k 



o 



' + !), fci 



(2^-l)(2£ + 3) 



(3.10) 



Alternatively, one can derive the coefficients k n using the monodromy technique explained in ap- 
pendix [Cj See appendix lC.il for the explicit computations. There are also asymptotic expansions 
of Kg for large aw; we refer the curious reader to [15] for further details. 

3.2 Singular Points 

We now focus our attention on global properties of equation (|3.7p . The branch points of solutions 
of the scalar wave equation for Kerr are located at r = r±,oo. The singularities at the horizons, 
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r±, are regular singular points while the singularity at infinity is an irregular singular point of rank 
1. One can identify the nature of the singular points, for example, by following the steps in section 
12.11 As it turns out in this case, all singularities in the ODE (|3.7|) are also branch points of the 
solutions. Another equivalent way of determining the nature of the singularity at, say, r = r+, 
would be to substitute a series expansion of the form 

R(r) = (r- r + ) ±ia + [l + 0(r - r+)] (3.11) 

and solve for a + and the coefficients of 0{r — r+). In this case, then, r + corresponds to a regular 
singularity. A similar computation for the inner horizon at r = r_ sets the value of the exponent 

a_ . We find that 

2Mr±uj-am 

a± := . (3.12) 

r + — r_ 

So we now know the conjugacy classes of the monodromy matrices associated to the horizons: 

/ e -27r« + q \ / e -27ra_ q \ 

M +"( o *•*)< M -^{ *--)■ (313> 

For irregular singular points of rank R, e.g., around r = oo, one must also include exponential 
factors in the series expansion. For our case, where the singularity has rank 1 at r = oo, the 
asymptotic expansion of the solutions to (|3.7[) is of the form 

R(r) = e =F*^ r =F<A-l ^ + o( r -i)] ; (3.14) 

where e ±27rA will correspond to the eigenvalues of the formal monodromy matrix, rather than to 
the eigenvalues of the true monodromy at infinity, e =t27rairr , as we will now explain in detail. 

At the singularity r — > oo, there will be a nontrivial conjugacy class of the monodromy matrix 
of the form 

( e ~2wa irr q \ 

^ J > ( 3 - 15 ) 

where a.i rr is yet to be determined and will be more challenging to compute than for regular singular 
points. Since infinity is a (non-resonant) irregular singularity of rank 1, there are also two Stokes 
matrices that must be determined, which we denote by Sk- Following the discussion in section [2^1 
this implies that the monodromy at infinity can be expressed as 

Moo = e 2niAo S^S . (3.16) 

The formal monodromy can be read off directly by replacing (|3.14p in (|3.7p . and recalling that 
circling r = oo in a positive direction means taking r — > e~ 2m r, 

( e~ 27rA 
e 2mAo =* ' ;. I , \ = 2Mu . (3.17) 
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0.15 



0.10 - 



<)M> 




Figure 1: The figure depicts the monodromy around the irregular singular point, ai rr , as a function 
of the frequency oj for fixed values of the mass M = 0.7, angular momenta a = 0.2, and t = m = 2. 
The black line outlines the analytic perturbative results for aj Tr while the grey line represents a fit 
to the numerical data given by the grey dots. 

It happens that A = a+ — a_ in the Kerr background, but we will only utilize this fact in our 
final expressions so that intermediate expressions are applicable to more general confluent Heun 
equations. In a basis where e 2m ° is diagonal, we define 



S-i 




,5'n 




(3.18) 



where Cq-i are the Stokes multipliers. Taking the trace of (|3.16p , we find that 



Oiir 



2tt 



cosh 



-i 



cosh(2vrA) + e^ A CoC_i/2 



(3.19) 



Computing oti rr is then equivalent to determining the product CqC-\, which is a rather involved 
task. Analytically, it is not obvious how to estimate aj rr for arbitrary values of the frequency. 
However, in appendix Owe develop a series expansion to compute ai rr directly for low frequencies. 
Alternatively, Stokes multipliers can be computed numerically to high accuracy using the method 
summarized in appendix [Bl developed by [16]. We developed the StokesNotebook [T7] (for Mathe- 
matica) that implements the aforementioned method to compute the Stokes multipliers and, hence, 
(|3.19p numerically. A comparison of the analytical and numerical results to order 0(w 4 ) in Kg and 
ai rr are shown in Fig. [TJ Naturally, the perturbative approximation breaks down for larger values 
of to. 



3.3 Scattering Coefficients 

For the purposes of this subsection, until we discuss quasinormal modes, the parameters in the 
radial wave equation (|3.7p are restricted to their typical physical ranges, i.e., < r_ < r+, w 6 R, 
t = 0, 1, . . ., m = —£, —£ + !,...,£, and hence Kg £ R. When we discuss quasinormal modes later, 
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the frequency u will be extended to complex values. 

Having identified the singular points, we are now ready to discuss reflection and transmission 
coefficients. As will be explicit below, global data (including Stokes data) is not sufficient to 
determine these physical observables. The derivations here are intended to quantify which data are 
fixed by the monodromy information, and which are not. We will return to this computation in 
section |U 

We have three singular points — inner and outer horizons, and infinity — hence from equation 
(|2.8p we know that 

M_M + Moo = 1 . (3.20) 

As we showed in section 12.31 finding the relation between two bases of solutions defined by their 
behaviors at regular singular points is largely algebraic. Still, we should specify clearly the choice 
of boundary conditions and its relation to the discussion in section 12. 3[ Therefore, our first task 
is to define ingoing and outgoing modes at the outer horizon and infinity, keeping the role of the 
monodromy as explicit as possible throughout the process. 

As in any other scattering problem, the phase of the solution should have a definite sign as 
the boundary is approached to correspond to waves falling into or out of the horizon. The ingo- 
ing/outgoing solutions at horizons are often written in so-called tortoise coordinates, r*, where they 
take the form of plane waves: 

i/> ~ e wr * (...) + e~ iuir * (...) , (3.21) 

(note that even though the solutions have essential singularities at the horizon, r* = — oo, we need 
not worry about Stokes phenomenon since the ODE, expressed in tortoise coordinates, will not 
be meromorphic in r* and the machinery we have discussed is not directly applicable in tortoise 
coordinates). For the outer horizon, adopting the left-eigenvectors of M + as a basis coincidences 
with ingoing/outgoing boundary conditions (|3.2ip . as the expansion in (|3.1ip indicates. The fun- 
damental matrix is therefore of the form 

$+ = 




for sign(u;a + ) = + , 

v + { (3.22) 

for sign(wa-| x 

where a+ is given by (j3.12[) . <!>+ are constant matrices, and <&^ is invertible. (When the parameters 
are instead complex, the sign of Re(w)Re(a + ) is the defining feature of ingoing versus outgoing.) 

At infinity, the plane wave basis of the form ()3.2ip diagonalizes the formal monodromy e 27r ° . 
Hence the basis at infinity, valid in an open wedge of the complex r plane that emanates from 
r = oo and contains the ray r > 0, is of the form 

"Ppw = ( Vin,pw ^out, pw J 
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(3.23) 



Here A = 2Mw as in (|3.17p . <J>n W are constant matrices, 3>g W is invertible, and ^ n r _n <J>™ w generally 
does not converge. 

Next, we can use the Wronskian (12. 4D to normalize solutions by denning an inner product on 
solutions via 

<^,y, 2 ) := W{ f^ 2) = ^'^'^ (Pi(r)arih(r)-Mr)Mi(rj) . (3.24) 

This is a constant because ipf will solve the same ODE as ip\ when the differential operator is real, 
as we have assumed for this subsection. We then define the norm of a solution to be ||^|| 2 := (ip,ip), 
which is proportional to the flux through a surface of constant r. In terms of this inner product, 
and under our assumptions on ranges of parameters (e.g., r + — r_ > 0), we find that 

, 9N + for 'Out' states, 

sign | M| 2 = (3.25) 

— tor In states. 

Translating this definition into the language of monodromies gives the following. When a mon- 
odromy a is real, the basis diagonalizing the monodromy matrix can be normalized so that the 
inner product defined by the Wronskian|£| which determines the flux through a fixed r surface via 
the Klein-Gordon inner product, is given by <t 3 |j This is the case in the Kerr background for a± 
at real frequencies, and it applies as well to the formal monodromy A. 

The computation of greybody factors (i.e., scattering coefficients) boils down to finding the 
connection matrix that transforms the two bases of solutions into one another. For our choice of 
boundary conditions at the horizon f|3.22|) and infinity (|3.23p . the connection matrix is defined as 

M pw ^+ = $ p >+ . (3.26) 

Since our basis can always be made orthonormal for q + and A real, the connection matrix can be 
brought to the form 

■M P w->+ = ( | : f J , \TZ\ 2 + \T\ 2 = 1 . (3.27) 

Hence, A4p W _>+ is an SU(1, 1) matrix, where \T\ and \1Z\ are the transmission and reflection 
coefficient, respectively. To compute (|3.2T|) . we will take advantage of the product relation (I3.20p . 
This product involves the true monodromy M m and not the formal monodromy, so we will split 
the computation into two parts: first we will connect $+ to the fundamental matrix $oo that 
diagonalizes M m via (|3.20p . and then we will relate $oo to $ pw using (|3.16p . 



As noted around (|2.4[) . recall that the Wronskian equals the determinant of <fr. 
6 Throughout, we will use a' to denote the Pauli matrices. 



16 



Note that the fundamental matrix that diagonalizes M^ can be described in Floquet form as 

(3.28) 




- oo 





In comparison with (|3.22|) . the expansion here is given by a Laurent series since r = oo is an 
irregular singularity. We emphasize that this implies that cti rr cannot be read off directly from the 
ODE. 

To convert from the basis of left-eigenvectors of M + to the plane-wave basis of left-eigenvectors 
of e 2n ° , we go through the intermediate basis diagonalizing M^ . Following section 12.31 the change 
of basis from M + to M m is given by 

/ sinh7r(aj rr -a + + a_) sinh7r(a irr + a + + a_) \ ,, nn , 

\ sinh7r(aj rr + a+ — a_) sinh7r(aj rT — a+ — a_) / 

where '~' denotes equal up to the equivalence relation (J2.15J) . This relates the two bases via 

$+ = $00^00^+ . (3.30) 

Through (|3.16p . the Stokes matrices tell us how to relate the basis of solutions diagonalizing the 
monodromy matrix M m to the plane wave basis that diagonalizes the formal monodromy e 2mA o, 
The Stokes matrices are given in (J3.18P and the formal monodromy in (|3.16p . From there, we can 
compute a change of basis between the two to be 

$oo = $p„.Mp„^oc ' ( 3 - 31 ) 

with 

sinh-7r(A — aj rr ) sinli7r(A + ai rr ) I 



(3.32) 



up to equivalence (|2.15p . 
From here we have 

M pw ^+ = M pw ^ 00 M 00 ^ + , (3.33) 

with the caveat that we have to determine the normalization of the solutions to make the product of 
two connection matrices meaningful. In particular, we have to determine the normalization of the 
columns in ^oo following (|3.24p . To start, for real frequencies, e 2nairr can be either real or a phased 
When e 27TOlirr is real, then the freedom in the normalizations of solutions can be partially fixed by 
making _A/f pw _i.oo and A^oo-s-+ into SU (1,1) matrices. When it is a phase, the normalizations of 



7 When the ODE is real, if ip is a solution then the complex conjugate ip* must be a solution as well. Furthermore, 
if rp has a definite monodromy e 27Ta ' rr then ip* has monodromy e 2jrc " irr . Since the two eigenvalues of M^ are e irr , 
it must be the case that e 2na ' rr is real or a phase. 
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the columns of ^oo can be chosen so that the inner product will be given by a . In this case, we 

Ipw- 



can require -M pw _>oo and .Moo->+ to be determinant 1 matrices that convert a cj 3 norm into a a 2 



norm, or vice versa. In fact, requiring the matrices to be SU(1, 1) for real e 27Tairr will automatically 
result in the latter property when it is a phase. 

Normalizing the columns of <&oo in this way allows us to write 





I 1 


Hi \ 




i! \ 


^ 


^ ) 


/da 


o d- 1 J 


■Rl 


1 


V ° d 2~ 




V 71 


Ti / 





AW = I 
where we have defined 



1 


no - 


\ / , 


—? 


TT 


1 ^ d^ 1 


■R.2 


1 


■h 


h / 


1 



Til = We" 



-2-7TQ,' 



sinh7r(A — ai„) 
sinh7r(A + airr) 



/sinh7r(a + + a_ + a, rr ) sinh7r(oq_ — a_ + aj rr ) 
y sinh7r(a + + a_ — aj rr ) sinh7r(a+ — a_ — aj rr ) 

7; 2 = 1-^ 2 for i = l,2. 

From the reality of the ODE — i.e., when the parameters such as u are real, as we've assumed for 
this section — we deduce that the parameters d\ and d% are phases while d,2 is a phase (respectively, 
real) when e 27rcw is real (respectively, a phase). 

With this parameterization, we can now read off the entries of .M pw _>.+ as reflection and trans- 
mission coefficients. The norm of the transmission coefficient is independent of d\ and d% in (|3.34|) . 
However it does depend on the unknown function cfe, adding another layer of complexity on top of 
the computation of the true monodromy a{ rr . Explicitly, we have 

\r\ 2 = i - \n\ 2 = — ^^2 r • ( 3 - 36 ) 

(d 2 + KiK 2 ) (d^ 2 + KiK 2 ) 

In summary, we cast the transmission coefficient as a function of Hip, which are specified in terms 
of the monodromies, and of the normalization d 2 . With the techniques presented here we cannot 
determine d 2 (or for that matter d\^). In principle this requires explicitly solving the ODE, which 
we have avoided. In section U we will discuss further how to constraint analytic properties of d 2 , 
while still avoiding solving the ODE. 

The advantage of this approach is that it highlights the dependence of the scattering coefficients 
on the global properties of the solutions. For one, it shows how the inner horizon data enters: even 
though we're studying a boundary- value problem on r G (r + , oo), the monodromy data of the inner 
horizon enters due to the global constraint (j2.8|l . This approach also appears to separate out the 
contributions that are intrinsic to the black hole from those that arise from the propagation to 
asymptotic infinity, and there is a sense in which this is true [9]. Although, as we have seen, the 
undetermined parameter d 2 plays a central role in understanding the transmission coefficient, and 
presumably mixes the intrinsic data of the black hole with that of asymptotic infinity. In section 
HI we will further constrain the dependence of d 2 on u), finding qualitative agreement with the 
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low- frequency approximations. 

3.4 Quasinormal Modes 

Quasinormal modes (QNMs) represent resonances in black hole scattering problems and therefore 
yield important information about the spectrum of radiation from a black hole. They are defined 
as solutions that are purely ingoing at the horizon r = r + and purely outgoing at infinity. For real 
frequencies, it is not possible to satisfy these boundary conditions, so we now relax that condition 
and allow u € C. This also means that the undetermined constants in (|3.34p . d\, d 2 , and d 3 , will 
become complex numbers, rather than just phases. 
In terms of the definitions in section 13.31 we have 

*in,+ = *out, P w , (3.37) 

which, along with the condition det.M pw ^+ = 1, implies that the connection matrix for a quasi- 
normal mode takes the form 

M«™ + = ( t T ) = ( l ~i ) . TV 0, (3.38) 

when Re(u;)Re(a+) < (the columns should be interchanged for the other sign). Since the fre- 
quency is complex, T' and 1Z' are generally not complex conjugates of T and 1Z. The quasinormal 
mode condition (I3.37p . or (13.381) . is viewed as a constraint on the frequency u; we denote those 
frequencies that obey (I3.38J) as w QNM . The expression in equation (I3.36J) is still valid, but with T' 
and TZ' replacing T* and TZ* , so TT' diverges at w = w QNM . 
From (|3.34p . we find explicitly 

_j_(d 1 d 3 (d 2 + d 2 ~ 1 iz 1 iz 2 ) ^(d 2 n 2 + d 2 - 1 n 1 ) \ 
PW " + W^^i + rt) jfc(<hKiR* + tf) ) • ( ^ 9j 

We want to bring this connection matrix to the form (|3.38p . The first constraint is that the product 
of the off-diagonal entries is equal to —1. Since the product is independent of d\ and c^, this can 
be rearranged to yield the equation 

(d| + 1Z x 1Z 2 ){d^ 2 + IZxlZi) = . (3.40) 



7ll^2 



Similarly, the product of the diagonal entries is independent of d\ and d% and must vanish (this is 
just the inverse of (J3.36P ). This implies 

■=^=i{4 + lZ{lZ 2 ){d^ 2 + 1Z{JZ 2 ) = . (3.41) 

Together with (|3.40p . this condition implies that {T\T 2 ) 2 / {1Z\1Z 2 ) either vanishes less rapidly than 
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the left-hand side of (|3.40|) . or does not vanish at all. 

Employing our results for the computation of the transmission coefficients (|3.35p and the data 
of the Wq N m from [2], we can easily see that T1T2 and IZ1IZ2 are finite and nonzero. Typical values 
of these constants can be found using StokesNotebook [17]; for Schwarzschild and Kerr black holes, 
see Tables [1] and [2] in appendix lB.il Appealing to the finiteness of 1Z\ and 7£-2 allows us to rewrite 
the condition (|3.40p more succinctly as 

dl + Tl{R,2 = or d 2 ~ 2 + TZiTZ 2 = , (3.42) 

naturally corresponding to the poles of the transmission coefficient (|3.35p . This computation makes 
it transparent that the spectrum of the black hole depends not only on the intrinsic data of the 
black hole, e.g., 7£i and IZ2 which are determined by the monodromies, but also on knowledge of 
the solution along r G (r+,00), which is encoded in d.2- 

That the frequency u) should satisfy at least one of the relations in (|3.42p is a necessary, but 
not sufficient, condition for that frequency to correspond to a quasinormal mode. The full set of 
solutions to (|3.42p correspond to those values of uj satisfying either in-out or out-in boundary con- 
ditions. To understand the appropriate subset corresponding only to in-out boundary conditions, 
we need to know the behavior of (dicfe) since its vanishing or divergence at solutions to (|3.42p can 
change whether the upper-left or lower-right entry of the scattering matrix vanishes|f| Without 
knowing efe, it is difficult to see this explicitly. In the example below for the interior resonances, 
the structure of the branches of solutions will be more clear. 

3.4.1 Resonances in the Interior 

A completely different boundary- value problem would be to impose boundary conditions at the 
inner and outer horizons. This setup has not attracted much attention, nevertheless it could be 
interesting for the purpose of studying the interior of the black hole, and since it is a boundary 
problem that we can solve precisely in our setup, it is worth discussing. 

Consider the connection matrix relating the left-eigenvectors of the monodromy matrices M + 
and M_ . Following our previous discussion (see section I2.3P , we have 

$ = § + M+^- , (3.43) 



8 A similar phenomenon happens for massive scalar fields in a BTZ background, which is solved by hypergeometric 
functions and can therefore be studied explicitly. The analog of (|3.42p would tell us that we have two branches of 
solutions indexed by integers ni,ri2 £ Z, and the restriction to in-out boundary conditions would then restrict m 
and n.2 to be natural numbers. If we wrote the scattering matrix in that case as 

di \ / ?f? ttt \ / d 3 



1 


R • 


T" 


T-> 


K; 


1 


T 2 


Tl 



^A! ijlo d .^ 

with IZ2 and 72 given by (|3.35[) , then we would find that (dicfe) would have poles or zeros at the solutions to the 
analog of (p£121) . 
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with 



,'ck o \ ( sinh7r(a + - a_ + a irr ) sinh-7r(a + + a_ + a irr ) \ / d 5 o 

• M +->- = ( o dj 1 J I . i / , v • i / s I I o rfr 1 ) » l 3 - 44 ) 

4 ' \ smn7r(a_|_ + a_ — aj rr ) smh7r(a + — a_ — aj 



l irrJ 



where d± and ofe are the undetermined normalizations of each basis. 

Consider imposing quasinormal- mode-like conditions on M-\ >_, given by (13.380 : i.e., we have 

\P; nj _ = >t out)+ . Requiring the matrix to be of the form (13.38P gives 

a+ — a_ + ai rr = in\ , rai £ Z , 

a + — a_ — a^r = iri2 , n 2 G Z . (3.45) 

We see, too, that the normalizations c^ and CZ5 contain information about the resonances in the 
connection matrix since we must tune them to match the form in (|3.38p . Depending on whether 
we choose in-out or out-in boundary conditions, there are two different branches of modes. The 
boundary conditions in-out and out-in are translated to regularity conditions of the solutions at 
each asymptotic region. Let's assume that Re(ia + ) > and Re(ia_) < 0, then the in-out boundary 
conditions will restrict the integers n\ t 2 in (|3.45p to be positive, and the out-in boundary conditions 
will correspond to n\ 2 being negative. 

These boundary conditions are analogous to those imposed on the angular equation (j3.6|) : see 
also appendix lC.il In the literature, these modes are known as spheroidal wave functions |18j . 

4 Symmetries of the Connection Matrix 

Having analyzed in detail the global properties of the ODE, we now turn to analyzing certain 
symmetries of these equations. By exploiting these very simple symmetries of the ODE, we will 
determine certain universal properties of the connection matrix Mi^j. We will mostly be inter- 
ested in finding Aii^j between a regular singular point and an irregular singular point, with the 
immediate applicability being to deduce further properties of the entries in ()3.34p . As we will show, 
this rather simple analysis has consequences for the spectrum of the black holes. 

To illustrate the construction, we first write explicit relations for the general class of confluent 
Heun equations (CHE) and then apply these to the problem we are interested in, namely the Kerr 
wave equations. The analysis can certainly be applied to other classes of ODEs, as in the case of 
the double confluent Heun equation (see, e.g., [18]) which is relevant for the Teukolsky equations 
in extremal Kerr. 

4.1 Symmetries of the Confluent Heun Equation 

This section follows |19| . using a more convenient notation for our purposes and correcting some 
typos. 

Consider the confluent Heun equation, with regular singular points at z = and z = 1, and an 
irregular singular point at z = 00. The regular singular points have monodromy eigenvalues e ±27ra ° 
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and e ±27rcn , respectively. The formal monodromy eigenvalues at infinity will be called e ±27rA , and 
the structure of the essential singularities at infinity will be given by e ±Jro2 . In terms of these 
parameters, the confluent Heun differential operator can be written as 



£[ao,ai,A,w,"](2) = Z ( Z ~ X ) 



d 2 f\-2ia Q l-2iai \ d /A . 



(ao + a\) 2 + i{pcQ + u\) + k iw(2iao + iX — 1) iw{2ia\ + iX — 1) 



z(z — 1) z z — 1 

The parameter k is defined such that when w = — i.e., when the singularity at infinity is regular 



then _ e ± 7ri v 1 +4>« are the monodromy eigenvalues. All parameters are generically complex|£| With 
the choice of scaling in the equation, the asymptotic behavior of solutions to £,\ aoaii x miK ](z)-ip(z) = 
are 

rp(z) ~ z 2iao , or 1 (z -> 0) , 

i(>(z) ~ (1 - z) 2iai , or 1 (z -> 1) , (4.2) 

V>(*) ~ e -Mtu* z iao4tei-iA-l ; ^ z ia +ia 1+ i\-l (^ ^ qq) . 

We could rescale VK- 2 ') = z* a °(l — z)* ai e~* ro;!; V( z )) i n which case ip(z) would have more symmetric 
asymptotic expansions around the singular points, but we keep the current scaling to more simply 
connect with the notation in [19]. When w = 0, the asymptotic expansion of solutions is convergent, 
and the leading behavior of an asymptotic expansion of the solutions around z = oo is 



iao +ia\ ± A \/T+4k— A, 



(4.3) 



On the other hand, when w ^ 0, the asymptotic expansion of solutions around z = oo is generically 
not convergent and k does not appear at leading order in the asymptotic expansion (not surprisingly, 
since the operations of expansion around z = oo and taking the limit w — > need not commute) . 
Let us first list the elementary symmetries of the confluent Heun equation, which relate the 
connection coefficients of two different confluent Heun equations to each other: 

£[ao,a 1 ,X,w,K}(z)-{z 2ia °'4>(z))=0 =► £[-a ,a 1 ,A, W ,«](*) ' fa) = ° > ( 4 - 4 ) 

A«o,ai,A,-A](^)-(( 1 -^ 2iai V ; (^))= => £[a ,- ai ,X,n,K)(z)-fa)=0, (4-5) 

-C[ao, ai ,A,«;,K] W • (e" 2iro ^(z)) = => £[ao, ai ,-A,- W ,K] W " fe = . (4.6) 



Denote a fundamental matrix that diagonalizes the monodromy matrix at z = z% by $ n m) I and 



,(*)0 



9 In relation to the Kerr radial equation (|3.7[) . we have that r = (r+ — r_)z + r_, ao = Q-, Qi = a+, A = 2Muj, 
zu = ixj(r+ — r_), and k — Ke + (a 2 — 8M )ui . Kerr is not the most general CHE since A = on — ao and since the 
parameters ranges are restricted. 

10 Note that the choice of <!>p W defined in (|2.29|) and section [3731 does not diagonalize the monodromy matrix Moo, 
rather it diagonalizes the formal monodromy (|3.16l) . Hence <3> pw diagonalizes e 2lTl °, while $o, 3>i, and <I>oo, diagonalize 
Mo, Mi, and M x , respectively. 



22 



the connection matrix (|2.14p between any two pairs of such fundamental solutions by M Zi -> z .. The 
elementary symmetries (|4.4[) and (|4.5p imply that the connection matrix between two regular points 
must be of the form 11 ! 

M = $ _1 $ =( fc ( a o,ai,A,w,«) fc(-a ,O!i,A,tJ7,/c) \ 

\ /c(ao) ~~ a i) A,n7, k) k(— (Xq, — a±, A, zu, k) I 



The same function k determines all entries of the connection matrix. Invariance under (|4.6p of <J>o 
and $1 implies that 

k(ao,a\, X,zu, k) = k(ao,a±, — X, — zu, k) . (4-8) 

Next, consider _A/f pw _>i, which relates the fundamental matrix diagonalizing Mi to the formal 
fundamental matrix diagonalizing the formal monodromy around the irregular singular point z = 
oo, (|3.16p . In a similar fashion, we can use (|4.5p and (|4.6|) to constrain the entries of the matrix: 

, / q(a ,ai,X,w,n) q(a , -«i, X,zv,k) \ 

M pw ^! = <P pw <Pi = . (4.9) 

^ q{a ,ai,-X,-zu,K) q(a ,-ax,-X,-zu,K) J 

The transformation (|4.4p further implies that 

g(aoi oil A, ro, k) = g(— ao,a\, X,zu, k) . (4-10) 

Another elementary symmetry arises from reflection z = 1 — u, which switches the two regular 
singular points. Call the generator of the reflection r. Let ip(u) := ip(l — u), then 

£[<*o,<*i,A, ro , K ](l -u) -V(l -u) = => ^[ ail Q ,A- ro ,«](") ■i'iu) = . (4.11) 

The reflection symmetry r has a particularly interesting implication. Consider how reflection acts 
on the infinite cover of the z-plane around the branch point at z = oo by writing v = 2z -i = P^ d > 
where £~R. Then z — > 1 — z sends v — > —v, which means 9^0 + (2n + l)ir for some neZ, This 
Z2 reflection, then, apparently lifts to an action of Z on the infinite-sheeted cover of the z-plane. 
In that case, we should consider the "reflection" operator r to act as r : (z, zu) — >■ (1 + e m z, e~ l7T zu) 
on the cover of the z-plane. 

Since $0 an d $1 have convergent series expansions, we can act with r term- by-term on their 
series expansions and deduce that 

r ' ^OjacctiAra.fslV 2 ) = $l,[ai,O!0,A,— w,k]\ z ) ) 

r • ^0,[a ,ai,A,ro,K]( 2; ) = ^Ojao^iA^^) > ( 4 -12) 

r ' ^ > l,[a ,ai,A,ro,K]( 2; ) = ^l,[a ,a!i,A,CT,/«]Oz) • 



One can also choose to normalize the fundamental matrices such that Mi^,j will have unit determinant. 
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Thus, acting with r on (|4.7|) implies that the connection matrix must satisfy 

Mx-+o(ao,a 1 ,\,zn,K)Mx-+o(ax,a ,\,-zn,K) = 1 2X 2 • (4.13) 

On the other hand, the action on $ pw cannot be naively deduced from the asymptotic expansion 
since it is not convergent and since the operation r 2 takes us to a different sheet of the z-plane. 
Instead, we see that 

r ' < ^pw,[ao,a 1 ,\,m,K,]{ z ) = < ^pw,[ao,a 1 ,\,e- 2,Ti 'Dj,K]\ e z ) 



-pw, 



B 



where in the last line we have used the fact that <I>p W and M^ are analytic in w[ 

In that case, then, there are interesting implications for the connection matrix between a regular 
and irregular singular point. Acting with r 2 on 

^X,la ,a u X,w,K]( z ) = ^pw,[a ,a u X,zu,n](z) M pw ^l(a ,ai,X,ZU,K) , (4.15) 

implies that 

M pw ->i(a ,ai, A, e" m w,K) = M OOt i a0tOiltX ,zu,K]M pw ^i(a ,ai, X,iu,k) , (4.16) 

where M^ is expressed in the plane- wave basis. So we see that the connection matrix has a branch 
cut in the complex ro-plane around w = with the same monodromy as <3? pw has around z = oo. 
Writing M m as 

Moo = e 2wiN ™ , (4.17) 

(see equation (|2.12p ). then a solution to (|4.16j) is 

M pw ^i =w~ N ^M(a ,a 1 ,X,w,K), (4.18) 

where M is a meromorphic function of w. The connection matrix is not a single-valued function 
of vj, and the branching structure is dictated by the monodromy at infinity via N^. 

The branch cut in the ro-plane that we point out here occurs in other contexts, as well. For 
example, it is also known to occur in the double confluent Heun equation [TB] (relevant to extremal 
Kerr), and even in the case of confluent hyper geometric (relevant to the quantum mechanics of 
a charged particle in a Coulomb potential) — in the latter case, it will not manifest itself in the 
transmission coefficient TT* for reasons we point out in section [4~2l 

By exploiting these very simple symmetries of the CHE, we were able to extract universal 



12 Since <!>p W is uniquely denned by its asymptotic expansion in a particular wedge of the z-plane, and since the 
coefficients of this asymptotic expansion will be analytic in vj, $ pw is analytic in -oo. By the definition of the 
monodromy matrix, analyticity of $ pw with respect to vj implies analyticity of M x . 
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properties of the connection matrix. As we will discuss in the following sections, these properties 
have implications on the spectrum of the black hole. In addition, there are also various interesting 
integral relations in [TS] that can be exploited to further constrain the entries of Ai Zi ^ z , though 
we will not use them here. 

4.2 Kerr Revisited 

As an application of the above discussion, let's see how these results constrain our expressions 
for the greybody factors and quasinormal modes of the Kerr wave equation. In section 13.31 we 
constructed the connection matrix, which in (|3.34p we expressed as 

AW+ = ( dl ^ ) -^pw^oo ( d2 d -i ) M^ + ( &A ,-x ) , (4.19) 

with A4 00 _ s . + and M. pvi -^oo given by (|3.29[) and ()3.32|) . and where d{ are undetermined parameters. 
On the other hand, the analytic property of the scattering matrix that we just deduced in (J4.18P 
states that 

AV->+ = ™- N °°M(a-,a + ,\,m,K e ) , (4.20) 

where w = uj(r + — r_). Since A^pw-^oo and A / I o->+ are meromorphic functions of w, the branch 
cut structure is necessarily encoded in the normalizations di in (|4.19p . To see this, first note that 



1 d- 1 ) -^Pw^oc (4-21) 

diagonalizes M m = e 2mNaa from the plane- wave basis and, hence, diagonalizes N^, too. So we can 
use (I4TT91) and (|4~20jl to write 

M(a-,a+,X,w,K t ) = ( dl ,-A M pw ^ («*"" _ J ( da d -A M^ + ( d;i d -A . (4.22) 



Since Ai(a-, q+, \,w,Kg) is meromorphic in w, this implies that 

d 2 = vj- ia ^d(a-,a+, A, w, K e ) , (4.23) 

where d(a_,a+, A,ro,iQ) is meromorphic in w. 

The implications of the symmetry transformations on the transmission coefficient are straight 
forward. From (|3.36|) we have 

so we see that the scattering coefficients are clearly not analytic in oj. In particular, a low-frequency 
expansion of the non-analyticity of (J4.23P gives 



*— j ^lOf-irr 



(u;(r + - r.)f [l - (2Muf ln{u(r + - r.)) ( pjgggpfcij ) + 0(lo* In a;)] , (4.25) 



25 



where we used the expansion of aj rr in (|C.36|) . The leading power of ui is in accord with results 
at low frequencies [20(, [2il I22|. with the caveat that we made a choice of branch in evaluating a, rr 
since cosh(27raj rr ) is really the well-defined quantity (so aj Tr is determined up to an overall sign 
and a shift by an imaginary integer). This ambiguity can be fixed in several ways, the simplest 
being by requiring the decay rate to be finite in the limit that Muj — > 0J 13 l 

A similar non-analyticity in the transmission coefficients has been reported in various places, 
for example, in [231 EH E3 El [8] . In particular, the authors in 0[8] constructed Floquet solutions at 
infinity (see (]3.28p ) and used them to match to the series solutions around the outer horizon. The 
challenge in using Floquet solutions is twofold: one, it is difficult to compute the Laurent coefficients 
without knowing a.{ rr exactly and without having enough boundary conditions for the coefficients, 
and two, the relationship to plane waves must be established to solve the desired boundary-value 
problem. The advantage, though, is that the Laurent series is convergent. The expressions the 
authors of [8] arrived at for the greybody factors are non-analytic functions of uj with the same 
branching structure we arrived at in f|4.23j) (|4.25|) from our simple analysis of symmetries of the 
CHE. 

Recall that the connection coefficient d 2 was particularly important in determining the location 
of quasinormal modes. With this new information, we find that ()3.42p becomes 

w" 2iCiirT d{a-,a + ,X,vj,Ktf +n 1 K 2 = or w 2iairr d(a-,a+, X,vj,K £ )- 2 +KiK 2 = . 

(4.26) 
Our simple analysis shows that the QNM spectrum of the Kerr black hole is governed by a tran- 
scendental equation. There are a few possibilities for solutions to this equation: 

1. Both terms separately vanish in (|4.26p . in which case if u is a QNM, then so is e 2m u. 

2. When u = Wq NM is a QNM, airr(ojQ NM ) is rational and both terms cancel each other. Then, 
the QNMs live on a finite-sheeted cover of the w-plane. 

3. Quasinormal modes require the cancellation of the two terms against each other and each 
QNM lives on a single sheet of an infinite-sheeted cover of the w-plane (oiirr(LUQ NM ) ^ Q). 

We can explicitly test whether the first two options are viable. It is clear from the data in Tables 
[1] and [2] that on rr is generically complex, and TL\R, 2 does not vanish for Wq NM . The final option is 
the conclusion we are inevitably led to. 

Further knowledge of how d(a_, a + , A, w, K() in (|4.26[) depends on its parameters is desirable. 
This would lead to a better physical understanding of the properties of the spectrum, such as the 
spacing of the modes, it could lead to new symmetries in the spectrum, or it could provide a simple 
way to see that all QNMs decay in time (i.e., that Im(w QNM ) < 0). 



The relationship between the greybody factor |7~| and the decay rate F is simply 

\T\ 2 2Mr+Lu - am 

1 = — ; , cx-\- = . 

e 47r Q+ _ i r+ -r- 
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The branch cut structure in the connection matrix can be attributed to the irregular singularity 
at infinity, making it difficult to disentangle which contributions of the QNM spectrum are intrinsic 
to the black hole. In the low-frequency limit, we see that there is a logarithmic contribution to 
the greybody that can be completely attributed to the presence of an irregular singularity [26|l27], 
Logarithmic branch cuts in the greybody factor have also been noticed, for instance, in certain 
zero-temperature correlation functions |24[ [28] . Since an extremal horizon generates an irregular 
singularity, this also would generate a logarithmic contribution to the greybody factor, fitting nicely 
with this zero-temperature result. This rather universal behavior in the transmission coefficients 
for extremal black holes has also been noticed in the analysis for specific examples of asymptotically 
AdS solutions, see e.g. [29] . 

We should emphasize that this sort of phenomenon is not particularly unusual in scattering 
theory. In fact, there is a familiar setup where these statements can be checked explicitly. Consider 
a scattering process governed by the confluent hyper geometric equation. We can then repeat the 
analysis for this case, except now the connection matrix .Moo-s-i- will just be the identity matrix 
since there are only two singular points. Still, A / I pw _>+ is non-trivial since infinity exhibits Stokes 
phenomenon. Then the non-analyticity in (|4.19j) arising from cfo corresponds to multiplication from 
the right by a diagonal matrix and, therefore, will cancel out when computing the transmission 
coefficient |T| 2 . Nevertheless, one can directly verify the branch cuts in the individual elements 
of the connection matrix by explicitly computing the scattering matrix and seeing that they are, 
indeed, present] 14 ! 

5 Future Directions 

In this paper we have developed a procedure, based around monodromy data, for understanding 
scattering coefficients and quasinormal modes in various black hole backgrounds. This development 
allowed us to demonstrate an important result concerning the analytic properties of the scattering 
coefficients in a Kerr background, also implying exact properties of the equations that govern quasi- 
normal frequencies. In arriving at this result, we utilized elementary symmetry transformations 
of the confluent Heun equation. However the ODE has additional integral symmetries which will 
further restrict the form of the connection matrix |314 [T5] ; it would be interesting to exploit these 
relations in more detail. 

It is rather challenging to obtain analytic properties of the quasinormal mode spectrum. Most 
of the work on quasinormal modes of rotating non-extremal black holes involves numerical methods 
[321 [33], though there are a few analytic methods for the Schwarzschild solution (e.g., [Ml ESI ES] ) 
and the (near-) extremal Kerr black hole (e.g., see |251 [37] and references within). Our results, 
valid in full generality, provide new insights into the complex-analytic properties of these modes. 



14 A simple example is to study scattering off a Coulomb potential, where this non-analyticity can again be verified. 
Finally, it's worth noting that the phenomenon of complex eigenvalues of a Hamiltonian, corresponding to bound 
states, being confined to a single sheet of the complex energy plane is a very generic feature of quantum mechanical 
problems; see, e.g., [30) . 
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A more detailed comparison with prior work, especially that on extremal black holes, would be 
intriguing. 

In fact, in some sense extremal Kerr could be more tractable than non-extremal: the two 
horizons merge and form an irregular singular point, yielding a double confluent Heun equation. 
The double confluent case has a much larger elementary symmetry group than the confluent case, 
and an analysis of its implications on scattering data has been nicely summarized in [18]. Appealing 
to results from [18] for a scalar probe, one learns that the scattering matrix has two branch cuts: 
one at u = 0, as for Kerr, and another at cj = ^tj, which is the critical frequency for superradiant 
scattering and is known to correspond to an accumulation point of quasinormal modes in the 
extremal limit [38]. We hope to return to a more detailed discussion of extremal Kerr in future 
work. 

We would like to emphazise that the analytic properties of the ODE governing the perturbations 
for higher spin modes is in the same universality class as the scalar probe, being described by 
confluent Heun equations; see appendix [El The technique used for the scalar field in section [3] 
will then apply for higher spin probes such as gravitational and electromagnetic perturbations. 
However, the problem of computing scattering coefficients has to be revisited. In particular, the 
boundary conditions and reality conditions of the modes might be modified. We leave these issues 
for future work. 

Of course, there are many other types of black holes that could be studied using these techniques, 
including charged black holes, those in AdS or dS spacetimes, and those in higher dimensions. In 
many cases, additional singular points appear at complex radii making the analysis more subtle, 
but in others, such as in AdS, asymptotic infinity is typically regular, so one does not have to 
worry about Stokes phenomena. While our exploration of scattering data in black holes focused 
on neutral, asymptotically flat backgrounds, we hope that the present paper will stimulate further 
progress in a much wider array of black hole backgrounds. 
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A Basic Notation 

Throughout this paper, we discuss ODEs of the form 




-O-UtHsW*- 



The notation we use to characterize different facets of the system, such as its solutions and global 
properties, are: 



• Mi'. Monodromy matrix, defined in (|2.6p . 

• e 2n °: Formal monodromy matrix, relevant for irregular singular points; see below (|2.24p . 

• $j: Fundamental matrix of solutions to (jA.ip . as defined in (|2.2j) . The subscript i could 
denote that the fundamental matrix diagonalizes Mj , as in (|2. 13[) , or it could denote one that 
diagonalizes the formal monodromy e 27r ° at an irregular singular point, as in (|2.23p . 

• Mi^j\ Connection matrix between <3?j and <F/, i.e., $~ $j, as in ()2.14p . 

• af Related to the eigenvalues of the monodromy matrix Mi, which are given by e \ Since 
tr(Mj) = 2cosh(2-7raj) is the well-defined quantity, cq is only defined up to an overall sign 
and a shift by an imaginary integer. 

• Aj: Similar to cti, related to the eigenvalues of the formal monodromy e 2 "" ° around an 
irregular singular point; see either (J2.23P or (IB.3p . 

B Estimating Stokes Multipliers 

To compute the Stokes multipliers numerically, we use the techniques developed in [16], which we 
summarize here. In the following, we focus on the case of a rank 1 singularity rather than a higher 
rank singularity, which can be found in |16j . 
Consider the differential equation 

d 2 Mz) + f(z)d x il>{z) + g(z)i>(z) = , (B.l) 

where f(z) and g(z) have convergent series expansions around z = oo£j We will assume that 
z = oo is an irregular singular point, so the solutions will have an essential singularity at z = oo. 
When the rank of the singularity is 1, we can expand f{z) and g(z) as 

OO j, oo 

/(*) = ££■ »(*) = ££> ( B - 2 ) 



15 In the case that f(z) and g(z) are rational functions, if we call h(z) the least common multiple of their denom- 
inators, it is computationally more efficient to multiply (|B.f [) by h(z) since it leads to simpler finite-order recursion 
relations than otherwise obtained. This is implemented in the Mathematica notebook we developed [T7] . 
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where rank equal to 1 implies that at least one of the coefficients /o, go, or g\, is nonzero. 

For simplicity we restrict to the non-resonant case where /q 7^ 4go, then the two formal series 
solutions to (IB.ip are given by 

00 00 

m*) = ^ Ai E ^ > fc(*) = eC2 ^ A2 E ^¥ > (B.3) 



s=o z s=o z 



where 



The exponents Ai 2 are given by 



Ci,2 = -^±M-go- (B.4) 



A i = —, -r- > A 2 = -7 ^- • (B.5) 

C2 - Ci si — C2 

The coefficients of the series expansion (IB.3J) are determined by the recurrence relations 

s 

(Ci -C2)so S) i = (s- Ai)(s- 1 - Ai)o s _i,i + ^ [Ci/j+i + ffj+i - (s-j- Xi)fj] ds-j-,1 , 

i=i 

s 

(C2 - Ci)sa s ,2 = - A 2 )(s - 1 - A 2 )a s _i i i + E [C2/./+1 + ffj+i - (s - j - X 2 )fj] a s - jy2 , (B.6) 

i=i 

with the starting values of ao,i = «o,2 = 1 and a Sj j = for s < 0. Generally, ^ s a s ^z~ s will have 
zero radius of convergence and only corresponds to an asymptotic expansion of a true solution 
around z = 00. 

In fact, to specify the asymptotic expansion of a true solution around z = 00, we have to 
specify the argument of z when we take the limit z — > 00 to determine what linear combination 
of formal solutions ipi and ip2 we will obtain. Thus, we can think of the formal solutions in (|B.3[) 
as corresponding to given true solutions only within a given wedge of the z-plane S k , with k 6 Z, 
defined by 

S k = {z : (k - i) vr - Arg (C2 - Ci) < arg(z) < (fc + ±) vr - Arg (C2 - Ci)} • (B.7) 

On a closed sector properly interior to S k ~i U 5fc U <Sfc+i, there is a single linear relation among the 
three associated solutions: 

ipk+x(z) = C k ip k (z) + Tpk-i(z) , (B.8) 

where the set of constants {C k }, with k E Z, are called Stokes multipliers (at infinity) of the 
differential equation (jB.ip . and the ip k f° r k ^ {1,2} are defined in (jB.lOp . Defining A k := 

e -(-l) fc fe(A 2 -A 1 )^ Cfc; we have that 

A k = { A ° f ° r k GVen ' (B.9) 

I A 1 fork odd. 
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We also have the relation 



1pk+2j(i 



Combining these, we see that 



3 2j\2Ki^ k ( e -2J7Ti z ^ 



for k even, 
for k odd. 



( Me~ 2ni z) Me~ 2m z) ) = ( M*) Mz) ) M °o > 



(B.10) 



(B.ll) 



where 



M n 



-2Ai7ri 





-2A27TI 



1 
C_i 



1 Co 
1 



(B.12) 



This is just as in (|3,16p . 

We can normalize the independent variable z in such a way that £2 — £i = 1 by replacing z with 

(B.13) 



C2-C1 



We implement this in the numerical approach. The values of the Stokes parameters will be related 
to the original by (£ 2 — Ci) ( Aa ~ Al ) times their old value. With the choice (|B.13|) . a good estimate 
of the Stokes multipliers is 



An 



-2iria 



8,2 



m— 1 

^ Oj,l r ( s - A 2 + Ai - j) 

j=0 



+ 0( S " m ), 



(B.14) 



A x = 2vri(-l) s " i a 



.s.l 



m—1 



Y, (-i)S-, 2 r( s + A 2 -Ai-i) 

i=o 



+ 0(s" 



For numerical implementation, this is a systematic way to compute A^i. Alternatively, if the 
coefficients a Sj i and a Si2 are known as a function of s, the Stokes multipliers can be obtained from 
the limit 



A ( 



2vri a s2 
hm — - 

00,1 s - > °° T(s — A 2 + Ai) 
2-7ri .. , _ „_-, a s 1 



\s-l 



A 1 = lim(-l) 

a s^oo 1 (s + A 2 — Ai / 



(B.15) 



When using the summation formula (|B.14p , which is a very efficient way to compute the Stokes 
multipliers, it's worth noting that the large s behavior of the coefficients a s ^ implied by (JB.15P 
suggests a restriction on the choice of m in (IB.14J) not mentioned in [16] . The sum over j in (1B.14|) 
would clearly diverge for m = 00, so one should truncate the sum around the smallest term. If we 
use the asymptotic behavior of a s i to get a rough idea where to truncate, this suggests choosing 



m 



< 



s + 1 



^Re(A 2 -Ai) 



(B.16) 
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n 


£ 


W QNM 


Txh 


n x n 2 


a irr 








0.220910 - 0.209791 i 


0.725516-1. 3456 i 


0.981555 + 0.183844 i 


0.10688-0.00489413?: 


1 





0.172234- 0.696105 i 


1.99733- 0.327733 i 


0.703811+0.213504?. 


0.184325- 0.395024 i 


2 





0.151484- 1.20016 i 


1.90982-0. 0988799 i 


0.638641+ 0.171458 i 


0.196496 + 0.139823i 



Table 1: This table shows the QNM of the Schwarzschild black hole (a = 0) of mass M 
n is the overtone and £ is the total angular momentum quantum number. 



2 where 



a/M 


i = m 


, ,Kerr 
W CMM 


K, 


7173 


K1K2 


Otirr 


0.2 





0.221535 - 0.209025; 


-0.0000179424 + 0.00030871; 


0.720005 - 1.34114; 


0.979579 + 0.17628; 


0.104914-0.0124874; 


0.4 





0.223398 - 0.206506; 


-0.0000966287 + 0.00123024; 


0.703135 - 1.32729; 


0.97394 + 0.152931; 


0.106939 - 0.0019500; 


0.6 





0.226342 - 0.201397; 


-0.000319103 + 0.00273531; 


0.673668 - 1.3024; 


0.966199 + 0.11148; 


0.106438 + 0.0022231; 


0.8 





0.229074 - 0.191402; 


-0.000841989 + 0.0046779; 


0.628351 - 1.2633; 


0.964207 + 0.0471877; 


0.10359 + 0.00897525; 


0.96 





0.222904 - 0.178774i 


-0.001356610 + 0.0061231; 


0.576185 - 1.22289; 


0.997748 - 0.00361162; 


0.094805 + 0.0129196; 


0.98 





0.221232 - 0.178961; 


-0.00134879 + 0.00633965; 


0.571018 - 1.21868; 


1.00052 - 0.00027528; 


0.0940499 + 0.0121101; 


0.99 





0.220894 - 0.178998; 


-0.00136306 + 0.00646113; 


0.568935 - 1.21602; 


0.999991 + 0.0000324433; 


0.0938947 + 0.0119539; 


0.999 





0.220768 - 0.178795; 


-0.00138927 + 0.00656798; 


0.569172 - 1.2139; 


1. + 2.44249 x 10~ lb ; 


0.0937403 + 0.0120083; 


0.9999 





0.218527 - 0.181398; 


-0.00123152 + 0.00660759; 


0.54043 - 1.21994; 


1. + 1.63758 x 10~ lb ; 


0.09376 + 0.00974293; 


0.9999 


2(11=1) 


1.98647 - 0.00705074; 


5.85717 + 0.00102783; 


1.00493 - 0.00035975H 


7.06779 x 10-" + 0.000264006i 


0.944334 + 0.498461; 


0.9999 


2(n=2) 


1.98644 - 0.021193; 


5.85719 + 0.00308939; 


1.00503 - 0.0107268; 


0.000207279 + 0.000347074; 


0.946068 + 0.495206; 


0.9999 


2(n=3) 


1.98635 - 0.0353136; 


5.85724 + 0.00514751; 


0.994687 - 0.0176253; 


0.000379994 + 0.000312636; 


0.949419 + 0.491461; 



Table 2: This table shows the QNM of the Kerr black hole, with mass M = ^, for different values 
of the black hole's angular momentum jj £ [0, 1] and values of £ = m. The overtone is n = 1 unless 
stated otherwise. 



where the upper sign is for Aq and the lower for A±. This is only a rough bound, so it is good to 
check different values of m and s to verify the consistency of any given computation. Obviously, 
this also implies a rough bound on the choice of s, namely 



s > max{l, 1 ± 2Re(A 2 - Ai)} 



(B.17) 



B.l Numerical Implementation for Schwarzschild and Kerr 

The above algorithm can be implemented in a straightforward manner for the ODEs relevant to 
Schwarzschild and Kerr backgrounds; the Mathematica code can be found in [IT]. For concreteness, 
we computed the monodromy data (|3.19p . (|3.35p . for certain values of quasinormal frequencies. In 
Tables [1] and [21 we used the values of w QNM reported in [2]. 

In addition to the low-lying quasinormal modes, we can also implemented the numerics for 
highly damped frequencies. Some interesting results on QNM for Schwarzschild black holes |34j 
showed that in the high overtone limit n — > oo, the frequencies are approximately given by Wq NM ~ 
Tb# m(3) + z7rTgfr(2n — 1), where Tbh = (87rM) _1 is the temperature. In this regime, we find 
that for large values of the eigenvalue £, the monodromy airr satisfies 



Re(a irr ) = aIm(w QNM ) + 6, a, b G 



(B.18) 



The plot in Fig. [2] illustrates this linear relation. 
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Figure 2: The figure shows the real part of the monodromy around the irregular singular point 
Re(aj rr ) as a function of the highly damped QNM frequency Im(w QNM ) = 7rTB#(2ra — 1) for fixed 
values of the mass M = ^i angular momenta a = 0, and £ = 1. The blue dots are the numerical 
values; for n 3> 40, the growth becomes linear and Re(ttj rr ) ~ 0.43243 Im(ojq NM ) — 5.51604. 

C Perturbative Monodromy at Irregular Singular Points 

Suppose we have an ODE 

(d z -A(z))$(z) = 0, (C.l) 

with an irregular singular point at z = oo, and suppose we can split the connection in two as 

4(z) = #(z) + e#(z), (C.2) 

where e is a small parameter and A^ ' (z) is a connection whose exact solutions are known, which we 
group into a fundamental matrix &°>(z). Then we can write &(z) = &(°> (z)^ 1 ' (z) , where ^^'(z) 
satisfies the ODE 

(d z - eA e )^{z) := (d z - e $(°)- 1 A( 1 )$(°))$( 1 )(z) = . (C.3) 

A t is what we would call the 'interaction-picture Hamiltonian' if we were talking about the Schrodinger 
equation (where A would be anti-Hermitian and <£' ) would be unitary). 

We want to compute the conjugacy class of the monodromy around z = oo, which is uniquely 
described by its trace and determinant. By construction, the determinant is one while the trace is 
given by 

trMoo = tr^z)- 1 ^- 2 ")) = tr^W^)^ 1 ^^^)^ 1 ^^^- 2 ")^^^ 2 ")) , (C.4) 

where we take z — > ze~ 2m to circle z = oo in the positive direction. Since we know exactly what 
&(°'(z) is, we can exactly compute the constant matrix 

Mg> := ^(zy^^ize- 2 ™) . (C.5) 
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Thus, we have 

trM^ = tr(MW$W( Ze - 27ri )$W(^)- 1 ) = ti(M^V{e e ^ Ae }) . (C.6) 



We can then expand the exponential and obtain trMoo as a series expansion in e. In fact, the 
right-hand side should be independent of a choice of 7 within a given homotopy class, so we can 
take 7 to be a circle centered on z = 00 and take the limit as z — > 00 to render the computations 
tractable. 

As we will see in both examples below, this method is very well suited for low-frequency expan- 
sions. It would be interesting to adapt a similar perturbative technique to high frequency regimes; 
in particular, it could be instructive to compare our technique with those used in [Ml ESI [39j 00] 
to compute highly damped quasinormal mode frequencies. 

C.l Modified Angular Equation: au <C 1 limit 

Introducing a coordinate z := cos 9, we can rewrite the angular equation (j3.6|) as 

d z {(l- z 2 )d z S(z)) + (K t + (acofz 2 - j^) S(z) , (C.7) 

which has an irregular singular point at z = 00 and two regular singular points at z = ±1 with 
monodromy eigenvalues ( — l) m . 

Before we begin computing the monodromy around z = 00, let's consider the other singularities 
first and understand the problem we are trying to solve. The monodromies around the regular 
singular points, z = ±1, are conjugacy equivalent to 

M±i = ( ( -J )m ( -Jr) • ( c - 8 ) 

Both solutions have a square-root branch cut when m is odd, but the nontrivial Jordan block 
implies more, it implies that one solution also has a logarithmic branch cut. We are looking for 
solutions that are regular at both z = 1 and z = —1, which correspond to 9 = 0, n. Since the 
regular solution at z = 1 (z = —1) corresponds to the only eigenvector of Mi (M_i), this means 
that Kg must take values that allow for M\ and M_i to share the same eigenvector. Recalling that 
we can write M±\ in a common basis as (see (J2.19P ) 

: „, ~ m I > M +1 =[ l> n ) , (C.9) 




*wr 



we see they share an eigenvector when 

e 2rra irr = j ==> a ^ g ^ and ivM ^ = 2 _ ( CJQ ) 

Computing the monodromy at z = 00 will then let us determine the values of Ki that satisfy this 
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condition. 

To compute M^, we will use (|C.6|) . We can write (|C.7p as a linear ODE with connection 

A(z) = ( ° 2 ^ )d* + (au;) 2 ( , ° M cfe 

=: ^ (0) (z) + (aw) 2 A (1) (z) , (C.ll) 

where we split K^ as 

K e = k + (au) 2 5K £ , k = £(£ + 1) . (C.12) 

The fundamental matrix for A^ ' can be written as 



Si(z) S 2 (z) 

(z 2 - l)d x S 1 {z) (z 2 - l)d z S 2 (z) 



^)=( , a ?£' ,, TJ„,, ) , (C13) 



where 



S 1 ( z ) = z &+m ~\(z 2 -\y^ 2 F 1 (a,^,l-5;z- 2 ) 

S 2 (z) = ^ +m "l(z 2 -irf 2 F 1 ( a + ^ + M + *;^ 2 ) (C.14) 



and 



In this case, then, 



5 = ^1 + Ak , a = i-f-i<5, p = h + a. (C.15) 



To obtain trMoo, we can now expand the right-hand side of (|C.6|) as a series in e := (auj) 2 so that 
toMx, =tr M£? + (au;) 2 i dztv{M^} A e {z)} + 0{{au:f) , (C.17) 

where A 6 (^) := ^(z)~ 1 A^ 1 \z)^(z). Calling trMoo = 2cosh(27ra irT ), this leads to 

2cosh(27raj r . r ) = —2 cos(ir yl + 4ko) 

-a 2 ^2vrsin(.yrT^) 3+4m2 ^-;^ 1+2 ^ ) + 0{(auf) . (C.18) 

Using ko = (■{£ +1), we have 

a irr = !! + ww^ — (Wi)A+i)-3) L + 0({auj)). (C.19) 

As we saw, to have a solution that is regular at both poles, we must set a.i rr € iZ, which 
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determines £ € Z and 



So, as we expect, Kg = ^n>o ^n( aw ) 2n where 

The restriction of — £ < m < £ comes from a more careful treatment of normalizability of the re- 
spective eigenfunctions than we gave here. One can continue this to higher order straightforwardly. 
For use in the next subsection, we quote the next coefficient from [15j . 



, _ (l-m-l)(l-m)(l+m-l)(l+m) _ (l-m+l)(l-m+2)(l+m+l)(l+m+2) /p 99 x 

K2 ~~ 2(2£-3)(2£-l) 3 (2^+l) 2(2m)(2<?+3) 3 (2£+5) ' l^-^J 

It is also worth mentioning that Mathematica has a built in function for computing .JQ(m, aw), 
called SpheroidalEigenvalue. 

C.2 Modified Radial Equation: Moo < 1 limit 

Recall the radial part of the scalar Kerr wave equation from (I3.7P : 



' AO (2Mr + w-am) 2 (2Mr„w - am) 2 , , „,. , rNN , 

(r — r_|_J(r_|_ — r_ J (r — r_)(r + — r_J 



fl(r) = K t R[r) . 



Next, define 

a := M\/l-e 2 =► r± = Af (1 ± e) , < e < 1 , (C.23) 

so e = 1 is Schwarzschild and e = is extremal Kerr (where the two horizons merge and become 
an irregular singular point). In terms of these parameters, the horizon monodromies are 



2Mr±oo - am 2Mw(l ± s) - m\f\ - e 2 .„ nAS 

a± := = . (G.24) 

r + — r_ 2e 



Next, for the non-extremal case, we definfi 

2Me 



+ r_, R(z) := z-^ 2 R(r(z)) , (C.25) 



16 For the computations in this section, we found it slightly more convenient to choose coordinates where the 
irregular singular point is at z — rather than at infinity. We then call its monodromy matrix M; rr to emphasize 
that it is the monodromy at the irregular singular point. This is what we called M x in other sections. 
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so that the equation becomes 

M1 - z)8 , + A. - i + jgl - ' + 4K ' - jj^g - 4e + «*> (C.26) 



4A/ 2 w 2 e 2 4eM 2 u; 2 (2 - el 



+ 3— 

z 6 



R(z) = , 



with the inner horizon at z = 00, the outer horizon at z = 1, and "asymptotic infinity" , the irregular 
singular point, at z = 0. We can turn this into a first order ODE with connection 

Mz) = ( 1+4.1 «j , ^-(^(7-4^) ^) +(2Mw)2 (£(^)-^ o)- (C - 27) 

Writing this as i(z) = i(°)(» + (2Mw) 2 iW(z), we choose A (0 \z) to be the first term in ([UT27]) . 
The solutions to the ODE with connection A^°>(z) are 

*<°>(z)=( ^ ** {z \ ), (C.28) 

V z(l-*W*i(*0 z(l-z)d z R 2 (z) J 

where 

Ai(z) = zl (\- z f ±S ^ 2 F x (a,P; 1 + 7; z) =: z* (1 - z)"^ 2 ^ F^z) , (C.29) 

« 7 a.-\-ft— 7 — 1 7 a + /3 — 7 — 1 

fl 2 (*) = *T*(1-*) 3 2 i ? i(a-7,/3-7;l- 7 ;z)=:z"2 (l _ z ) 2 F 2 (z) . 

This satisfies 

det$ (0) (z) = -7 (C.30) 

for all z. The constants appearing as arguments of the hypergeometric function are given by 

a = — y~ + i(a+ + a_) , /3 = —^- + i(a+ — a_) , 



7 = y/1 + 4FQ - (2Mw) 2 (7 - 4e + e 2 ) . (C.31) 

Similar to the angular equation, because the transformation R = \[zR has branch cuts, the trace 
of the monodromy for the original system (j3.7|) differs from that of A(z) by an overall minus sign. 
Thus, 

Mfl. = -m{z)^m{z^) = - ( ( " o 1)7 ^ ) = -e^ 3 . (C.32) 

Now we want to begin computing the Dyson series for the connection 

A e (z) := ^(z)-l ( _ £2 _%_*) ° ) *«»(*) , (C.33) 

where we identify the e expansion parameter (not to be confused with e) with (2Muj) 2 . Writing 
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this in a slightly more useful way, 



e(e + (2-e)z) _ 2ia+ ^ I -F 1 (z)F 2 (z) -F 2 (z) 2 \ ^ 
1^ \ F 1 (z) 2 F 1 (z)F 2 (z) 



We must keep in mind the validity of the order to which we work: each term in the Dyson series 
comes with a power of (2Mw) 2 . On the other hand, M and to appear in the solution & > as well, 
so the answer is only meaningful after expanding these constants to the appropriate power of w. 
Finally, 

tiM irr = -2costt7 + (2Mu;) 2 j> dz tr (M^Aq(z)) 

+\(2Muf ldz 2 idz 1 tv^M^.V{A n (z 2 )A n (z 1 )}^+0(u J 5 ) . (C.35) 

Keeping terms of order no more than w 4 , we find 

a 2 irT := (^cosh-^itrM^)) = -£ 2 + ^ a n (2Muj) n , (C.36) 

n=l 

where the first few coefficients for £ ^ are 

01 = , 

15^ + 1) -11 

a 2 = £- 



' (2£ + 3)(2l + l)(2l - 1) ' 

r - 7 51(1 + 1) -3 



+ l)(2* + 3)(2£ + l)(2*-l) ' 

« = mi +!)(«+ IK* + ff + vm - vm - 3) ( 80(f {l + 1))5 [115 + 2fe2] 

-8(1(1 + l)) 4 [6971 + 450^ + 1234e 2 - 1020m 2 (l - e 2 )] 

+ (£(£ + l)) 3 [97075 + 21480^ + 8058e 2 - 38928m 2 (l - e 2 )] 

-(1(1 + l)) 2 [358211 + 6(11827 + 60e 2 - 5529m 2 (l - e 2 ))] 

+3£(£ + 1) [7387 + 7854^ - 252e 2 - 702?n 2 (l - e 2 )] 

-45[49 + 121£ + 72m 2 (l-e 2 )]) . (C.37) 



The leading term when £ = m = is instead given by 

a 2 rr = -f (2Mu) A + 0(u 5 ) . (C.38) 

The expression for the monodromy is also valid in the limit e — )■ 1, which corresponds to the 
Schwarzschild solution, though the intermediate steps require a slight modification because the 
parameters of the hyper geometric functions a and /3 are equal, requiring special treatment. The 
limit e — > may require a separate treatment since the coordinate transformation (JC.24P is singular; 
although, the low-frequency expansion we obtain above certainly has a well-defined limit e — > 0, so 
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it may also give the correct expression for the monodromy in the extremal case. 

Another approach to compute a.{ rr was discussed in [7], where they reported only the leading 
correction in u). Their discussion involved studying the recursion relations for a Floquet solution 
of the type (13.28J) . A careful analysis of the low- frequency limit allows one to use their results to 
solve for a% r up to 0(w 2 ), which agree with the values of a\ and a 2 computed above. 

D Schwarzschild Black Hole 



In Schwarzschild coordinates we have 
ds 2 = 



1 + ?M.\ d f + ( 1 _ ^M_) l dr 2 + r 2 d0 2 + r 2 gin 2 M ,2 _ 



(D.l) 

The Klein-Gordon equation for a massless scalar is again fj3.4f) . and expanding in eigenmodes 

i/>(t, r, 6, 4>) = e- iu)t+im<t) R{r)S{6) , (D.2) 

the spherical function S(0) satisfies 

1 



sm( 



d 9 (sm8d e ) 



m 



2 i 



sin 



S(6) = -K e S(9) 



(D.3) 



whose solution is the associated Legendre polynomial P / m (cos 9) with eigenvalue K^ = 1(1 + 1). The 
radial equation for R(r) is 



d r r(r- 2M)d r + 7 ^rpr + r 2 + 2M(r + 2M w 2 

(r — 2M) 



R(r) = K e R(r) . 



(D.4) 



The singular points of the ODE for Schwarzschild corresponds to the coordinate singularity r = 0, 
the horizon r = 2M, and infinity r — > 00. 

To determine the scattering coefficients, we first must find the constants a+, a_, and aj rr , which 
are logarithms of the eigenvalues of the monodromy matrices associated with each of the singular 
points (see equations f)3. 13j) and (|3.15p ). We can find this information by studying the behavior of 
the solutions R(r) near these points. Around the regular singular points r = and r = 2M, the 
behavior of the solutions is respectively 



R{r) = (r - 2M) ia + (01 + ...) + (r - 2M)~ ia + (a 2 + . . .) 



Q- 



2Mu 



R{r) 



n+b 3 \ogr + •••)+ r~ ia -{b 2 + ---) , a_ =0 



(D.5) 



Note that the singularity at r = is a resonant regular singular point (hence the logarithm). The 
singularity at r — > 00 is an irregular singular point of rank 1, and the problem of finding a !rr will 
be analogous to that of Kerr, so it can be determined by (|3.19j) . While the product of the Stokes 
multipliers CqC-± will differ, the formal monodromy will be the same as that for Kerr, A = 2Muj. 



39 



E Higher Spin Perturbations 

Define eigenmodes for a spin-s field ip s as 

Mt,r,e,(j>) = e^ t+m ^R s (r)S s (9) , (E.l) 

which are related to scalar (s = 0), electromagnetic (s = ±1), and gravitational (s = ±2) pertur- 
bations around the Kerr background. The detailed relation between ip s and the field perturbations 
can be found in [20| l2T] . 

For s = 0, ±1,±2, we define R s (r) := A~ s ' 2 R s (r). Then the radial part of the wave equation 
for a particle of spin s can be written as 

* a a , i 2Mujr + ~ § s ( r + - r -) ~ am ) (2Mwr_ + §s(r + - r_) - am) 

a r /\a r -\ - [hj.z) 

[r — r + )(r + — r_) [r — r_j(r+ — r_j 



+uj 2 r 2 + 2(Mw + is)ujr + 2Mw(2Mw - is) - s 2 - K e 



R,(r) = . 



Again, K^ s is the separation constant, i.e., it corresponds to an eigenvalue of the spin- weighted 
spheroidal equation 



1 „ ,. .„, , „ , 2 (m + s cos 8Y 2 
-d e (sin Ode) + (acj cos 6 - s) 2 - — — - s 2 + K e 



sin ti sin 



S s {0) = • 



Notice that the radial ODE is not real when s/0; instead, complex conjugation (for real parameters 
and real r) is only a symmetry if coupled with (a, m, s) — > (—a, —m, —s). 

The form of the radial equation is the same as before: a confluent Heun equation with regular 
singular points at the two horizons and an irregular singular point of rank 1 at infinity, so the same 
techniques from the spin analysis carry over here. The monodromies of R around the horizons 

are given by 

is 2Muj — am ,_ , 

a± = =f- + . (E.3) 

2 r + — r_ 

Around infinity, 

so in the notation of section 14.11 we have 

w = oj(r + -rJ) , X = 2Muj + is, k = K i>s + s 2 + (a 2 - 8M 2 )w 2 . (E.5) 

Notice that a + ± a_ differ from their s = values by an imaginary integer, since s E Z, and 
similarly for A, so reality properties of quantities such as sinli7r(a + ± &_) are unaffected. 

Massless fermionic perturbations (s = ±1/2) need a separate treatment |4H 142], Furthermore, 
reality conditions and boundary conditions are sightly different in this case. However, an expression 
like (|R3|) still applies. 
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